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Abstract 

CO 

*^ ■ We reconsider, from a novel perspective, how unitarity constrains 

OO 

the corrections to the ratio of shear viscosity r\ to entropy density 
s. We start with higher-derivative extensions of Einstein gravity in 
asymptotically anti-de Sitter spacetimes. It is assumed that these 
S-h ' theories are derived from string theory and thus have a unitary UV 

completion that is dual to a unitary, UV-complete boundary gauge 
theory. We then propose that the gravitational perturbations about 
a solution of the UV complete theory are described by an effective 
theory whose linearized equations of motion have at most two time 
derivatives. Our proposal leads to a concrete prescription for the cal- 
culation of rj/s for theories of gravity with arbitrary higher-derivative 
corrections. The resulting ratio can take on values above or below 
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1/4-7T and is consistent with all the previous calculations, even though 
our reasoning is substantially different. For the purpose of calculat- 
ing rj/s, our proposal also leads to only two possible candidates for 
the effective two-derivative theory: Einstein and Gauss-Bonnet grav- 
ity. The distinction between the two is that Einstein gravity satisfies 
the equivalence principle, and so its graviton correlation functions 
are polarization independent, whereas the Gauss-Bonnet theory has 
polarization-dependent correlation functions. We discuss the gravi- 
ton three-point functions in this context and explain how these can 
provide additional information on the value of rj/s. 
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1 Introduction 



The gauge-gravity duality [JJ [2] can be used to learn about the hydrodynamic 
properties of strongly coupled gauge theories in the limit that their number 
of colors N is infinite. This is achieved through the study of weakly coupled 
anti-de Sitter (AdS) gravity theories of one dimension higher [31 H]. 

In the case that the theory in the AdS bulk is Einstein gravity, the ratio 
of the shear viscosity (rj) to the entropy density (s) in the dual gauge theory 
takes on the well-known value of [5j El [3j [7] 

2 - r- m 

s 4tt 

in natural units. Kovtun, Son and Starinets [8] conjectured that 2 > j- 
holds for all strongly coupled fluids, and this restriction has since been re- 
ferred to as the KSS bound. 

We are considering generic extensions of Einstein gravity in an (asymp- 
totically) AdS spacetime of dimensionality D > 5. Our focus is on static 
black brane backgrounds, relevant for studying the hydrodynamics of the 
field theory (FT) dual. We assume that there are higher-derivative terms in 
the extended theory, which amounts to finite iV or quantum corrections to 
the FT. The number of derivatives can be arbitrary, but the extensions are to 
be regarded as perturbatively small relative to the Einstein term. Our focus 
is on corrections to the ratio of two-point functions rj/s; however, as a conse- 
quence of our results, we find that quantities depending on the higher-point 
functions can similarly be studied. 

The entropy density s can be directly obtained from the Wald formula 
[9j [10] , whereas the viscosity 77 can be determined (at least in principle) by a 
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suitable generalization of the Kubo formalism. The calculation of the ratio 
rj/s for gravity theories having two-derivative linearized equations of motion 
is well understood. We refer to such theories as "two-derivative gravity" or 
sometimes as the "Lovelock class". These are the theories that belong to 
the special class of Lovelock theories [TT] or are, like f(7Z) gravity, closely 
related. 

But, if the linearized equations of motion contain terms with four or more 
derivatives, the theory inevitably contains Planck-scale ghosts and it is not 
clear whether these ghosts affect the value of the hydrodynamic parameters. 
So the situation is more complicated than often appreciated, and this matter 
has not yet been fully resolved. This is the aim of the current paper. 

This issue of ghosts has been discussed in [T2], where rj/s was calculated 
for the FT duals of four-derivative gravity theories. The ghost modes and 
their influence on the computed values of rj/s were explicitly exposed and 
understood. Then a set of boundary conditions (BC's) that eliminates the 
ghost modes was proposed. It was further claimed that, once these BC's are 
enforced, the validity of the KSS bound is ensured. 

The conclusion of [12J can be confronted with the particular case of 
Gauss-Bonnet gravity. This is a special theory for which the Lagrangian 
contains four-derivative terms but the linearized equations of motion contain 
at most two derivatives, and so it is two-derivative gravity by our definition. 
For this case, it is well known that the ratio rj/s can be below the KSS 
bound. Apparently, by having no ghosts to begin with, the Gauss-Bonnet 
model is able to evade the need for enforcing the BC's of [12] while violating 
the bound. 
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Here, we extend the considerations of [12] and propose that metric per- 
turbations about a solution of string theory are always described by two- 
derivative gravity This should apply in particular, to black branes in AdS 
space. So that, given a higher-derivative theory, the resulting effective two- 
derivative theory should be used for all calculations of the two-, three- (and 
higher-) point functions. This includes (but is not limited to) s and rj, as 
these parameters are directly related to the two-point functions of specific 
gravitons in the black brane background. 

This restriction to two-derivative gravity, for which there is a limited 
number of choices, implies that all such gravity theories and their FT duals 
depend on only a small number of parameters. Moreover, the above pro- 
posal already provides a clear distinction in philosophy between what we 
are suggesting and the previous methods for calculating r]/s. The resulting 
two- derivative theory should not be viewed as a perturbative extension of 
Einstein gravity (even if its higher-derivative origin is) since one uses the 
information from the UV-complete parent theory to allow the propagator to 
change. And so the arguments for enforcing the KSS bound in [12] no longer 
apply. 

The rest of the paper proceeds as follow. Section 2 starts with an expla- 
nation on how to extract the viscosity from the one-particle irreducible (1PI) 
effective action, with special attention to the implications of higher-order 
corrections. We then propose a novel set of BC's that maintains explicit uni- 
tarity by reducing any higher-derivative theory to an effective two-derivative 
model. It is made clear how our proposal differs from previous methods, 
which are either explicitly or implicitly based upon a Wilsonian form of the 
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effective action. Next, the entropy density is considered; in particular, the 
application of the standard calculations of s to a higher-derivative theory. 
We are, at the end, able to prescribe the ratio rj/s for a higher- derivative 
theory in an unambiguous way. 

In Section 3, we explain that, for current purposes, any such two- derivative 
model must be either an Einstein theory or a Gauss-Bonnet theory. The 
value of rj/s is shown to be sensitive to this choice and we provide a crite- 
rion to determine which theory should be used. This criterion relies on the 
distinction between the Einstein and Gauss-Bonnet forms of the two-, three- 
(and higher-) point functions. 

The final section briefly summarizes our conclusions and explains their 
relevance to multi-particle correlations in strongly coupled fluids. 

Appendix A contains the details of a calculation presented in Section 3, 
while Appendices B and C review relevant known results. 

2 A prescription for evaluating rj/s 

We are considering a theory of gravity with an Einstein-Hilbert term and 
arbitrary higher-derivative corrections. The main objectives of this section 
are to explain how the 1PI effective action can be used to extract the viscosity, 
propose BC's that ensure unitarity and provide a precise prescription for 
calculating rj/s. An important conclusion of this section is that any such 
gravitational theory should be a two-derivative gravity. Although our current 
focus is primarily on the viscosity (and, then, 77/s), a similar procedure can 
be applied to any graviton n-point function. 
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2.1 Initial considerations 

We assume here that the theory describing the AdS bulk spacetime is string 
theory, a unitary and UV-complete theory. In the low-energy limit, as rele- 
vant for calculating transport coefficients like rj, string theory can be approx- 
imated by an effective gravitational action, which includes the lowest-order 
(Einstein) term and possibly some higher- derivative corrections. 

We further assume that the full UV-complete theory has a black brane 
solution. We do not need to know the exact form of this solution, only its 
geometry. 

As a concrete example, we have in mind the case of a .D- dimensional AdS 
black brane for which the geometry is described by the following background 
metric 

ds 2 = -F{r)dt 2 + + ^dx 2 . (2) 

The index denotes the transverse space dimensions i — 1, . . . , D — 2, Lis 
the AdS radius of curvature and r is the radial coordinate (orthogonal to the 
brane). The function F vanishes on the horizon at r = r^, F(rh) = , and 
asymptotes to r 2 / L 2 as the AdS boundary (r — > oo) is approached. 

According to the gauge-gravity duality, the bulk has a UV-complete and 
unitary FT dual that "lives" at the AdS boundary; for example, M = 4 
supersymmetric Yang-Mills theory [1] or the Kats-Petrov class of M = 2 
conformal (C)FT's 

We wish to study small metric perturbations about the solution of the 
full theory and describe them in terms of an effective field theory. Since 
our interest is in physical observables in the bulk, we should calculate on- 
shell graviton scattering amplitudes near the boundary of the AdS spacetime. 
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These amplitudes are then dual to gauge-invariant correlation functions of the 
energy-momentum tensor of the FT, whereas other types of bulk amplitudes 
are generally not physical. 

Let us make a cautionary comment about the on-shell condition in AdS 
space. Since the physical quantities are really on-shell boundary amplitudes, 
it should have been possible to describe them purely in terms of boundary 
quantities. However, in practice (as discussed in detail below), one needs to 
define the solutions of the linearized equations also in terms of bulk quanti- 
ties. For instance, in a black brane background, one of the boundary condi- 
tions is usually imposed on the horizon. This involves gauge choices in the 
bulk, and one should make sure that the final results are physical and not 
affected by any of these choices. 

The relevant effective bulk action for the gravitons is, therefore, the 1PI 
action as derived from string theory, since this serves to generate the physical 
on-shell amplitudes. This 1PI action comes already equipped with an on-shell 
condition that defines the solution of the linearized equation for gravitons 
around a black brane background. 

2.2 Extracting 77 from the 1PI action 

Let us recall how the 1PI action is derived from string theory, starting with 
the case of M = 1 supergravity in 10-dimensional Minkowski spacetime. The 
on-shell condition for such theories is 

□e V = > (3) 
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up to quite high order (as explained shortly). D E is meant as the differential 
operator that satisfies the equations of motion from an Einstein theory. Cor- 
rections to the 1PI action start with terms that are quartic in the Riemann 
tensor (specifically, quartic contractions of the Weyl tensor [IH[T5]) and so 
contain 8-derivative terms. These terms determine only the 1PI four-point 
function (and possibly some higher-point functions) while not changing the 
one-, two- and three-point functions. This is the standard outcome when 
evaluating the 1PI action order by order in perturbation theory; new inter- 
actions do not change the lower-point functions. In particular, they do not 
change the one-point function (equivalently, the on-shell condition) nor the 
two-point function (equivalently, the propagator). 

When one considers the 1PI effective action in a lower-dimensional AdS 
Schwarzschild geometry (e.g., D = 5), the situation changes. Because the 
spacetime is no longer flat, the Weyl tensor is no longer vanishing in the new 
geometry. So that, if one computes the corrections coming from the quartic 
(Weyl 4 ) terms, the on-shell condition, propagator and three-point function 
are all modified. That the on-shell condition and the propagator are changed 
by interaction-induced corrections is a consequence of treating the full theory 
in a non-perturbative way rather than perturbatively near a specific vacuum 
solution. For instance, terms that only contribute to interactions in a given 
vacuum can modify the propagator in another. This is the main difference 
between the perturbative treatment discussed in [T2] and the current non- 
perturbative treatment. 

The FT viscosity rj is a physical quantity and is, therefore, defined by a 
physical on-shell amplitude of a specific class of perturbations, the transverse- 
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traceless gravitons h xy . Q| Of course, the viscosity rj can be denned, in prin- 
ciple, purely in terms of a FT prescription that is completely independent of 
the bulk. 

Let us now review how the on-shell condition is enforced in a black brane 
geometry for the case of the bulk theory being Einstein gravity. Then we 
will discuss the on-shell condition for more complicated situations. 

First, the Einstein equation for the gravitons is linearized in some gauge. 
The conventional choice is the radial gauge, which sets all graviton pertur- 
bations h^ r to zero such that is arbitrary. The resulting linearized form 
of the Einstein equation is DeIi^ = . 

The next step is to choose the BC's. The standard choice is Dirichlet 
BC's on the boundary of AdS and incoming BC's at the horizon of the black 
brane. These choices completely fix the solution of the linearized equations 
and, thereby, uniquely determine the on-shell conditions for gravitons in this 
background. 

The solutions of the linearized equations are of the plane-wave form 
h^u ~ h ilv {r)e lPX ~ %VlT , where X are the space coordinates and T is the 
time coordinate for both the FT and the brane. P is the dimensionless space 
momentum and Q is the dimensionless frequency; both of which are in units 
of the brane temperature. 

One then uses the solution of the linearized equation to evaluate the on- 
shell value of the two-point function of gravitons near the AdS boundary. By 

1 The labels x and y denote transverse brane directions that are mutually orthogonal 

as well as orthogonal to the direction of propagation. 

2 In the radial gauge, the {x, y} gravitons decouple from the other polarizations and, 

so, this gauge is compatible with transverse and traceless h xy 's. 
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the standard rules of the AdS / CFT duality, this correlator is related to the 
gauge-invariant two-point function of the energy-momentum tensor of the 
FT. Since the FT (as well as the bulk) is in a thermal state, this on-shell 
two-point function is also the 1PI retarded propagator. 

The viscosity rj and the other transport coefficients of the FT can be 
extracted from the two-point function of the energy-momentum tensor or, 
equivalently, from that of the gravitons near the AdS boundary in the hy- 
drodynamic limit. In this limit, the frequency Q and the momentum P both 
vanish, possibly at different rates. 

In practice, the value of rj can be evaluated by a number of means. For 
instance, via the Kubo formula [5] , the membrane paradigm j3] or simply by 
looking at the two-point function (h xy h xy ) at the horizon (equivalently, at 
the associated kinetic terms in either the action or the equations of motion) 
[T6l [T7J [TS]. The latter is so because, in the hydrodynamic limit and when 
the equations of motion contain at most two derivatives, the effective action 
has the same value at the horizon and at the boundary. 

2.3 Extracting rj from higher-derivative gravity 

Let us now discuss how the on-shell condition is enforced when the lowest- 
order Einstein-Hilbert action is corrected by higher-order terms (such as the 
previously discussed quartic terms). To fix the on-shell condition here, we 
will use exactly the same principles that were used above in the simpler case. 
A slight complication arises when the on-shell condition gets corrected at 
higher orders. In which case, there is some additional freedom in defining 
the on-shell condition due to field redefinitions adding terms proportional to 
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the lowest-order equations. This freedom will play a significant role in what 
follows. 

To see that the on-shell condition for the gravitons is modified by the 
presence of higher-order corrections, we can look at the linearized equation 
for the h xy gravitons. As discussed in [TU] and clarified in [TjJ], the La- 
grangian of a general extension of Einstein gravity can be expressed in terms 
of the metric, the Riemann tensor and its symmetrized covariant derivatives, 
while anti-symmetric combinations of derivatives can be traded off for more 
Riemann tensors, £ = C [g ab , 1Z abcd , V ai TZ abcd , V \ ai V ' a 2 )1Z a bcd, ■ ■ ■] ■ The lin- 
earized equation for the h xy gravitons can be obtained by expanding the full 
equations of motion 

dC 1 

- 2V a V b X" b + Tl abcp X ahc q - -g pq C = (4) 



to linear order in h xy . Here, 



dTZabcd 1 ai Tlabcd\ 

+ V (ai V a2) ( ® C J + ... (5) 

\C[V( ai V a2 ) ><<abcd\ J 

and the ellipsis means ever-increasing numbers of symmetrized derivatives. 

The linearized equation can also be regarded as a perturbative expansion 
in derivatives. Schematically, suppressing tensorial structures and indices, 

E + ec 2 (r)VV + eLc 3 (r)(V) 3 + eL 2 c 4 (r)(V) 4 + . . . 



+ e n - 1 c 2n _ 1 (r)L 2 "- 3 (V) 2 "- 1 + e"- 1 c 2n (r)L 2 "- 2 (V) 2 ™+ ...jh xy = 0, 

(6) 
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where O e is the operator for the Einstein equation ( n E h a b = ), e <C 1 is 
a dimensionless perturbative parameter controlling the strength of the cor- 
rections and the dots denote increasing numbers of derivatives as well as 
mass *™ of pe rtrub a tl ve orde, S Tie te^s, *<r), *(r), e te „ M e a 
model-dependent collection of dimensionless radial tensors that are at most 
of order unity but can have subleading contributions of order e (and higher) 
from yet higher-derivative corrections. These tensors represent the most 
general way to contract a given number of covariant derivatives and are 
built with the background metric, background Riemann tensor and its sym- 
metrized derivatives. Consistency of the perturbative expansion fixes the 
lowest power of e for a given number of derivatives. For classical a' string 
corrections, which correspond to finite iV corrections for the FT, e ~ lp/L 2 
(with l p being the Planck length). On the horizon, reduces to □, the 
standard d'Alambertian. 

Let us recall from [12] how this procedure works in the simplest case of 
four-derivative corrections to Einstein gravity. Eq. ([6]) then has the simpler 
form (up to henceforth implied mass terms) 

U E h xy + eb(r)D E h xy + eL 2 a(r)D 2 E h xy = 0. (7) 

To summarize the preceding discussion, the linearized equation for the 
{x, y} gravitons is generically of the form 

(n E + eV+...)h xy = , (8) 

3 One can absorb such mass terms into the definition of so that these are essentially 

irrelevant. Further, one can disregard perturbative terms with a single derivative acting 

on a graviton, as these are related to mass terms by integration by parts. 

4 This is only strictly true in the hydrodynamic limit. See |12] for clarification. 
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where T> denotes the leading-order higher- derivative corrections to the lin- 
earized equation (coming from the leading corrections to the Einstein-Hilbert 
action). The exact form of T> will not be essential for what follows. 

The operator T> generically contains terms with four or more derivatives, 
including some with at least four time derivatives. The equation will include 
more and more time derivatives as the perturbative order increases. So that, 
to define the on-shell condition, Eq. ([8]) has to be supplemented with new 
BC's in addition to those already specified in the Einstein case. Recall that 
some of the BC's are fixed at the boundary of AdS and others on the horizon 
of the black brane. The most relevant BC's for the current discussion are 
those at the outer boundary. 

That Eq. contains more than two time derivatives is very problem- 
atic because of the ghost modes that must appear in such cases. Since our 
treatment is classical, the problem manifests itself in severe instabilities if the 
ghost modes have non-vanishing amplitudes. (See [12] for additional discus- 
sion on this point.) What makes the situation a little more bearable is that 
the ghosts are typically at the cutoff (Planck) scale and, therefore, expected 
to be irrelevant to the calculation of hydrodynamic coefficients like rj. 

This situation when ghosts are present but have a cutoff-scale mass is 
ambiguous in connection to the calculation of rj (and any other coupling). 
One might reason that it is safe to disregard their presence. However, to 
proceed, one should in any case ensure that the ghosts have decoupled before 
the computation is carried out. Rather than rely on a general argument that 
may or may not apply, we propose to exorcize the ghosts from the get go; 
ensuring that they do not have the opportunity to interfere with the low- 
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energy physics. 

In fact, by the same line of reasoning, it is safer to eradicate all the 
extra modes with cutoff-scale masses, ghosts or otherwise. We propose to 
accomplish this goal through a particular choice of BC's. The rational for 
this choice is that, if some properties of the cutoff-scale modes determine a 
low-energy physical quantity, then the resulting value of this quantity cannot 
be trusted. So that, whenever a calculation involves the cutoff-scale modes, 
one should make sure that their amplitude, couplings, etc., do not enter into 
the final result. The best way to achieve this is, in our opinion, to set them 
to zero from the very start. 

We are now ready to make our central assumption about how the BC's 
should be chosen. We propose that 

The boundary conditions imposed on the higher- 
derivative equation have to be chosen in a unique way 
such that, from the many solutions to the higher-order 
equation only one mode has a non-vanishing ampli- 

(9) 

tude. This mode is the unique mode that, as e vanishes, 
connects continuously to the solution of the lowest-order 
equation (T3|) when supplemented with its original bound- 
ary conditions. 

The essential difference between the proposal (Q and the proposal in [12] 
is that, here, we are proposing to set to zero the amplitude of all the cutoff 
scale modes and not just the ghosts. For four-derivative theories, the two 
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proposals are equivalent. 

Let us now describe how the proposal (Q is implemented for gravitons 
in an AdS black brane background in the presence of corrections to Einstein 
gravity. To implement our choice of BC's on Eq. (JBJ), let us suppose (for 
simplicity) that, at leading order in e, the operator T> is obtained from a four- 
derivative or Riemann 2 correction. In [12], it was shown that it is possible 
to factor Eq. ([8]) in the following way: 

D E + tV = [1 + tb 2 {r) + ea(r)D E ] [1 + e&i(r)] D E , (10) 

where a(r) and b(r) = ^(r) + 6 2 (?") are functions of the AdS radial coordinate 
that are uniquely determined by the operator T>. To verify this form, one 
need only expand out the right-hand side to first order in e and compare it 
with Eq. ©.I 

From Eq. f lTU]) . we can extract the two modes that the corrected theory 
contains (now suppressing graviton indices): one massless graviton h\ and 
one additional mode hi whose mass is near the cutoff scale 1/e. Only one 
additional mode is added in this case because the highest derivative is of 
fourth order. The distinction between the modes is still subject to the afore- 
mentioned freedom of field redefinition. Here, the freedom is manifested in 
the functions b\ and &2 5 which are arbitrary provided that their sum remains 
fixed. This freedom cannot affect physical quantities and indeed it does not. 

For a fixed identification of modes, the BC's are imposed according to (J9]) 

5 Thcre are subtleties regarding the commutation of derivatives; see [12] . Alternatively, 
one can bypass this issue by imposing the limit of large (radial) momenta on the gravitons, 
so that the coefficients are effectively constant. 
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as follows: The BC's for the massive solution h 



2, 



[1 + eb 2 (r) + ea(r)D E ]h 2 = 



(11) 



are h 2 = at the boundary and incoming at the horizon. This constrains the 
amplitude of h 2 to vanish everywhere in the bulk (and obviously at the AdS 
boundary) by setting both its normalizable and non-normalizable modes to 
zero at the boundary. For the massless graviton, one chooses the "standard 
BC's", incoming at the brane horizon and Dirichlet at the AdS boundary. 

The 1PI effective action to a given order in gravitons is obtained by 
expanding the original higher- derivative effective action to the same order 
and, then, setting the extra mode to zero as required from the solution of its 
equations of motion. The result is a two- derivative action for the (modified) 
massless graviton. It is still not uniquely defined because of the freedom to 
change b±, while keeping b± + b 2 fixed. Ultimately, though, physical quantities 
must not depend on this freedom (as argued above). 

The situation in which the equation contains terms with higher than 
four derivatives is conceptually the same as we have discussed, although 
technically more complicated. The factorization procedure is expected to 
follow along similar lines and produce an additional factor for any additional 
pair of derivatives. Out of all the modes, one still has to choose the single 
mode that connects continuously to the massless graviton and set BC's such 
that all the rest vanish. The theory must then necessarily reduce to quadratic 
order or 



where the b(r) here is completely analogous (and equally ambiguous) to the 
coupling in the four-derivative case. We have not, however, implemented 



(1 + eb(r)) Dsh = , 



(12) 
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the proposal in a concrete example derived from a theory with six or more 
derivatives. 

Once the excess modes have been eradicated (and the freedom in b has 
been resolved), what is left is a two-derivative theory for which the usual 
AdS/CFT rules apply. One determines any n-point function by turning on 
a source (the non-normalizable mode) for the "normal" graviton, taking the 
n'th derivative of the action with respect to this source and then evaluating on 
the AdS boundary (with some prescription for removing divergences). Since 
77 is expressed in terms of a 1PI two-point function, which is also the two- 
point function of the FT energy-momentum tensor, the proposed procedure 
then computes 7] in an unambiguous way. 

The very same reasoning extends to all other transport coefficients and 
higher-point functions. 

The main conclusion of this section is that, when the proposed BC's are 
enforced, the resulting 1PI effective action for the remaining mode is a two- 
derivative gravity (meaning that the linearized equations of motion derived 
from it contain at most two derivatives). 

2.4 Extracting 77 from the Wilsonian effective action 

Another way of solving Eq. (JSJ) and calculating 77 is via the Wilsonian effective 
action. In practice, the calculations are only presented at the level of the 
equations of motion and done in the "reduced-action method" of Banerjee 
and Dutta [20J. All other accepted methods (e.g., [2T]-|25J) are equivalent 
to it. 

By this method, one iteratively reduces the higher-derivative equation to 
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a two-derivative form as follows: The lowest-order field equation 02]) can be 
written explicitly as (suppressing graviton indices) 

[d 2 + f 1 (r)d r + f 2 (r)]h(r) = 0. (13) 

The functions f\ and f 2 depend on the background and can also depend on 
P, Q and T. Their explicit form is not relevant for this discussion. 

Eq. (fT3|) is then used iteratively to reduce the order of the operator T> in 
Eq. OS]) until the full equation reduces to a quadratic equation. For example, 
when the operator T> has a four-derivative term as its higher order, then one 
has to use Eq. QTB"]) twice. Clearly, the procedure amounts to integrating out 
the massive modes, and so the equations that the massless graviton obeys 
can be derived from the corresponding Wilsonian effective action. 

In our notation, each iterative step is equivalent to using the relation 

n E h = 0{e)h. (14) 

Then, for example, 

D 2 E h = D E [D E h] n E [0{e)h) ~ eU E h. (15) 

One does this until there are at most two- derivative terms left and then 
discards any subleading term in e. The final result of this iterative process 
is then to reduce the higher-derivative Eq. (jSJ) to a quadratic in derivatives. 

For instance, when applied to the field equation <^7§ of a four-derivative 
corrected theory, this reduction process gives back simply [T2] 

(l + eb{r))n E h = . (16) 

Eq. QT6]) is then solved, as usual, with Dirichlet BC's at the AdS boundary 
and incoming BC's at the black brane horizon. Eq. ffTB"]) can also be utilized 
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to find a quadratic (Wilsonian) effective action, which may then be used 
(according to the standard AdS/CFT rules) to calculate bulk amplitudes 
near the boundary of AdS or, equivalently, FT correlation functions. 
In short, one ends up with [12] 

77 oc 1 + eb{r h ) , (17) 

which can be verified by extrapolating the usual two-derivative methods. 
(The constant of proportionality is fixed by Einstein's theory.) But, im- 
portantly, b(r) is the exact same b(r) = b\{r) + &2(?~) that appeared in the 
non-unique decomposition in Eq. (JTUJ). 

2.5 Comparison of the 1PI and Wilsonian approaches 

The Wilsonian calculation can be compared explicitly to the 1PI calculation 
with our proposed BC fl9]). The result is that the Wilsonian approach leads to 
a solution with a non-vanishing, order-e amplitude for the redundant mode. 

The two procedures give the same equation for the gravitons up to order 
e, and only differ at order e 2 . Meaning that, if one wishes to calculate to 
order e, it is possible to use either method. Since, in general, the Wilsonian 
effective action is not well defined in string theory, the justification for its use 
comes from the fact that it is expected to (and indeed does) lead to results 
that are equal to the ones obtained from the 1PI effective action at each 
order in perturbation theory. However, at higher orders this may become a 
delicate balancing act. 

The final result is that both methods agree on the value of t] (ambiguity 
aside) to leading order in e. This is somewhat of a red herring, as both 77 
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and s are densities and, as such, depend on the choice of coordinates, field 
definitions and so forth. Only the ratio rj/s is an invariant physical quantity. 
Meaning that it is now time to bring s into the discussion, which we do next. 



2.6 Calculating the entropy density 

The entropy density s is, similarly to rj, related to a graviton two-point 
function [2"oT, IT9"] . and so its calculation is subject to some of the same con- 
siderations. 

Let us momentarily consider a two-derivative theory of gravity that is not 
necessarily Einstein's theory. It is by now established CEO HE] that the 
extraction of rj amounts to reading off the horizon value of b in an equation 
formally identical to f fl6|) . Perhaps less well known is that one can play 
the exact same game to extract the entropy S or its density s. In this 
case, the gravitons in question are the h rt gravitons, which decouple from 
all others at the horizon^ The precise equivalence of this procedure to the 
Wald formulation was first pointed out in [26J and then made explicit in [19J. 
That is, one can look at the h rt component of the field equation 



and read off the horizon value b s (rh). This represents the correction to the 

entropy density when the Einstein value has been normalized to unity. 

6 The zero-modes of the h rt gravitons do not vanish on the horizon despite the choice of 
radial gauge. This gauge is inconsistent with the identity h a b — S7 a Xb + V&Xa (x a is the 
Killing vector), which has to hold on the horizon [9l[T0]. This is similar to the zero-mode 
problem in the Coulomb Gauge in electromagnetism and is resolved in a similar way. 
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A subtle difference in the 77 and s calculations becomes apparent only 
when higher- derivative corrections are introduced. Recall that 77 has an ar- 
bitrariness that is introduced by the presence of a ghost mode. On the other 
hand, s is somehow "smarter". As clarified in [19], any term with four or 
more derivatives can be reduced to two derivatives (or less) via successive 
applications of the Killing identity, V c V a X& = —T^abcdX d ■ This reduction 
process follows from a relation between the h r t gravitons and the horizon 
Killing vector \ a that is valid only at the horizon (see footnote E]) • 

It follows that, at the horizon r = r^, the {r, t} counterpart of Eq. ([6]) for 
a higher-derivative theory is 

{l + ec 2 {r h ))Uh rt + ... = , (19) 

where the additional terms that are represented by the ellipsis are only e- 
order mass terms. Higher-derivative terms do not exist! This simplification 
is unique to this class of perturbations and makes the entropy density special 
amongst graviton couplings. 

Strictly speaking, as we have just argued, the Wald entropy density s 
is insensitive to higher-order corrections and to the ghost modes that they 
induce. Thus, s appears to be free of the type of ambiguities that haunt r\ 
and other couplings. However, in the general context of the discussion, this 
conclusion is misleading because it holds only at the horizon, [j As we move 
away from the horizon, the ghost modes and their associated instabilities 
come back. 



7 Although the Wald density can only be denned at the horizon, one can use FT ther- 
modynamics and AdS/CFT rules to define s in terms of the h t t and h Xi x . gravitons. This 
definition has been shown to limit to the Wald formula as r — > Wf\ . 
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Our conclusion is that, to define rj or s or any other physical coupling, 
the only safe way to proceed is to restrict the discussion to two-derivative 
theories. 



2.7 Calculating the ratio r)/s 

We have tentatively identified the viscosity for an effective two-derivative 
theory as rj oc 1 + efe(r/ l ). However, only a ratio of densities, such as rj/s is 
physical. The ratio rj/s is gauge invariant and also invariant under bulk field 
redefinitions. The most convenient way to evaluate it is, as explained in [28] 
and reviewed in Appendix B, to use a particular choice of field redefinition. 
This choice calibrates s for a given theory to its Einstein value when compared 
at fixed values of or (equivalently) at fixed temperature. So that, in this 
field-redefinition frame, s is set to its Einstein value irrespective of the of the 
theory, whether corrected or uncorrected by higher-derivative terms. 
Using the proposed field redefinition, we find that 



s s E 



eb{r h ) = - 1 - l + eb{r h ) , (20) 



47T 

where a subscript E denotes the Einstein value and the tilde on b indicates 
that it is the previously defined coefficient but calculated in the field-redefined 
theory. 

Let us note that this leading-order result is in agreement with the previ- 
ous methods (as discussed in Subsections 2.4-2.5), although our reasoning is 
substantially different. 
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3 Two- derivative effective actions of higher- 
derivative gravity 



The focus of this section is on answering the following question: What is the 
two-derivative theory of gravity that effectively describes, by way of an action 
principle, a given higher- derivative model? The answer involves an important 
observation about the relationship between two-, three- (and higher-) point 
functions and how this relationship can be used to further constrain the ratio 

T]/S. 

3.1 Initial considerations 

To begin, let us recall the effective two- derivative field equation for the {x, y} 
gravitons, 

1 + efc(r)] D E h xy = . (21) 

We wish to find a covariant action from which this equation is derived. Before 
delving into the technical details, let us try to guess what the answer should 
look like. 

As argued earlier, the effective theory should be a two-derivative covari- 
ant theory of small gravitational perturbations about a background solution. 
This is already enough to tell us that it must be one of the Lovelock theories. 
Then, depending on the dimensionality of spacetime, there are only a few 
possibilities. For D = 5, the choices are just Einstein and Gauss-Bonnet 
gravity. These considerations apply to the two- and three-point functions 
and, in a limited sense to be clarified later on, also to the higher-point func- 
tions. Further simplification occurs when restricting attention to the calcula- 
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tion of rj/s because, for any dimensionality, Einstein and Gauss-Bonnet are 
then the only relevant theories. This was observed in [16] and is explained 
in Appendix C. Each of the two theories leads to a distinct value of rj/s. 



3.2 A tale of two theories 

Let us first consider the Einstein case. We wish to find an "Einstein-like" 
Lagrangian whose variation gives Eq. ( l2"Tj) . The answer is 



« inffi] . (22) 



Here, b(r) can be viewed as a radially dependent correction to the gravita- 
tional coupling or an energy-scale dependent coupling from the perspective 
of the dual FT. That is, 



— > 



1 + eb(r) , (23) 



16nG D 16irG D (r) 16%G D (r -> oo) 

where Go is Newton's constant and we have used that higher-derivative 
corrections die off when r — > oo . 

The actual sense in which the coupling in (I22p is radially dependent is 
as follows: First, the radius r is fixed to some specific value; for instance, 
the boundary of AdS or the horizon of the brane. Then one calculates the 
1PI amplitudes for that chosen radius. The amplitudes would then be in 
agreement to that of Einstein's theory with the standard fixed Newton's 
constant. 

In spite of appearances, the variation of Eq. ( 122|) leads to Eq. (T2Tj) without 
involving derivatives of b. Any resulting term in the variation of Eq. (12 2p 
contains two or less derivatives, and so a derivative on b means one or zero 
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derivatives on the graviton. Such terms are then (irrelevant) e-order mass 
terms. [] 

Now, one finds via the standard (two-derivative) methods that t] and s 
are both corrected by the same factor of 1 + ebfrh), and so their ratio is (of 
course) the usual Einstein value 



= ^ • (24) 
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For the Gauss-Bonnet case, the calculation is technically more involved. 
As explained in Appendix A, Eq. (T2Tj) would follow from the variation of 

L GB oc TZ + — 

+ eL 2 (3(r) [TZ 2 - ATZ ab TZ ab + TZ ahcd TZ abcd ] , (25) 

where /3 is a radial function such that, in the proximity of the brane horizon, 

P(r) = — ^ Mr) . (26) 

f K ' 2(D - l)(D-A) y 1 K ' 

Again, (3 can be viewed as a radially dependent coupling parameter that 
should be fixed at a radius of choice. 

One then finds that r] goes again as 1 + e6(r^). However, now s is uncor- 
rected, as it would be for any Lovelock theory with an unmodified Einstein- 
Hilbert term. (This follows from the argument of Appendix C, which is 
indifferent to a radially dependent coupling.) Hence, 

-1 = ^ [l + <*(r*)l • (27) 

Depending on the sign of eb(rh) , the ratio [tj/s]gb can be either smaller or 
larger than the Einstein value. 



8 Any term in the Lagrangian whose variation gives rise to a single derivative on a 
graviton is, up to surface contributions, equivalent to a mass term. 
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3.3 The two, three (and higher) point functions 

Since Einstein and Gauss-Bonnet gravity lead to different hydrodynamics, 
one may ask how else the two theories are distinguished and whether such a 
distinction can be further used to determine the ratio rj/s. 

Let us first suppose that there are no higher Lovelock terms, as would 
be the case for D = 5. (We comment on Lovelock extensions at the end.) 
Then, once the value of 77/s has been determined to a given order in e, 
the graviton three-point function is determined to the very same order (as 
clarified below). In particular, the n-point functions of Einstein gravity are 
polarization independent, as follows from the equivalence principle, whereas 
the n-point functions of a Gauss-Bonnet theory are polarization dependent 
in a specific way. Indeed, this is the very same polarization dependence that 
allows the value of rj/s to deviate from the Einstein value l/An. 

For the case of two-derivative theories, it was argued in [2SJ ED] that 
the two- and three-point functions are determined by a small number of 
interaction terms. Further, the two-point function is already enough to fix 
the three-point function, provided that the central charges of the FT dual 
are known [3T]. The premise behind these arguments can be understood as 
follows: The quota of at most two derivatives per term and general covariance 
conspire to severely restrict how one can modify the propagator. If mass 
terms are disregarded (these can anyhow be absorbed in a redefinition of 
□e), a modification of the linearized equations of motion is restricted to 
□ e — >• + J-" ab VaVb , where J-" a6 parametrizes the allowed corrections and 
its tensorial structure is determined by the background geometry. 

We conclude that this intricate relation between the propagator and three- 
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point function can be used to constrain the choice of two-derivative theory. 
A detailed presentation of the Einstein and Gauss-Bonnet two- and three- 
point functions (as well as their higher-point functions) will be deferred to a 
subsequent publication, as the exact expressions are not immediately useful. 
Their basic forms, however, can be deduced from general covariance, the 
two-derivative limit on the linearized field equations, and the applicability of 
field redefinitions and simplifying gauge conditions. 

The end result is that the Einstein and Gauss-Bonnet two-point functions 
can respectively be expressed in terms of a few coefficients, 

(h ac h bd ) E = Mfg ah g cd h ac h M + M E 2) g ab g cd g ef V e h ac V f h M , (28) 

(h ac h bd ) GB = (M^ B g ab 9 cd 9 ef + ^ 2) B n cedf g ab + vg B 1l acbd g ef ) V e h ac V f h bd 
+ (MP B g ab g cd + QP B U acbd ) h ac h bd . (29) 

Here, the coefficients M and Q are functions of r and the rest are num- 
bers. The superscript or 2 denotes the number of derivatives acting on the 
gravitons. 

Similarly, the three-point functions can be expressed as 

(h ac h be h df ) E = A^g ab g cd g ef h ac h be h df + g ab g cd g ef g^ h ac V l h be V 3 h df , 

(30) 

(h ac h be h df ) GB = {A^ B g ab g cd g ei g^ + B ( * B K^g ab g cd + C%\ll c ^g ab g^ h ac V t h be ^M 
+ (A^ B g ab 9 cd 9 ef + vf B n ced f g ab ) h ac h be h df , (31) 

where the coefficients A and V are functions of r and the rest are numbers. 



28 



To simplify, we have used that, for a theory of pure gravity, the Ricci 
tensor goes as [32J lZ a b = ^1ZS a b and, for Einstein gravity only, the Rie- 
mann tensor reduces to TZ ab cd = D ^_^ Tl [S a c 5 b d — S a d 5 b c ~\ . Further, we 
have dropped any terms that are redundant via symmetries and assumed the 
transverse-traceless gauge (extending to a general gauge would add more 
terms but is otherwise straightforward). It is also possible to express any 
n-point function for these particular theories using similar expressions, how- 
ever, we will not do so in this paper. 

The crucial point is that the Riemann tensor appears only in the Gauss- 
Bonnet correlators. The Einstein and Gauss-Bonnet n-point functions are, 
therefore, distinguishable in way that cannot be eliminated by a field redefi- 
nition. 

Now, given any higher-derivative theory, the simplest way to find the 
correct two-derivative description is to perform a field redefinition so as to 
bring the theory to the appropriate Lovelock one. When the choice is strictly 
between Einstein and Gauss-Bonnet, it becomes the simple matter of inspect- 
ing the higher- derivative action or field equation for polarization-dependent 
terms or, equivalently, for the appearance of the four-index Riemann ten- 
sor lZ abcd . If none appear, the theory has a natural completion to Einstein 
gravity and r]/s reverts to its Einstein value of 1/Att. Otherwise, the theory 
naturally completes to Gauss-Bonnet gravity, leading to the value for r]/s as 
in Eq. (j2Tj). 

Both of these outcomes are compatible with our previous prescription ( 120]) 
in Section 2. This agreement is trivial when the choice is Gauss-Bonnet, so 
let us consider Einstein gravity which is a polarization-independent theory. 
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Then, any field redefinition that changes s must change all two-point cou- 
plings by the same amount. Hence, for the Einstein case, s — > se implies 
that i] — > i]e or b — > b = 0. 

For the purpose of evaluating higher-point functions, additional input is 
necessary. This is on account of possible Lovelock extensions when D > 6 
and, irrespective of D, information about the higher-point functions that is 
lost in the linearized theory. So that, in general, further information besides 
polarization dependence versus independence is needed to clarify the precise 
identity of a polarization-dependent theory. Such additional input could be 
a perturbative hierarchy between the coefficients of the interaction terms. 
For example, suppose that the terms quadratic in Riemann tensor are of 
order e while those quartic in the Riemann tensor are of order e 2 . Then 
the magnitude of the deviation of the ratio r]/s from its Einstein value 1/47T 
and, similarly, the magnitude of the deviation of the three-point function 
from its Einstein value can be used to find out which term is responsible 
for the leading correction. Once this leading contributor is identified, then 
its contribution to four- and higher-point functions in some kinematic region 
can also be determined to leading order in e. 

4 Conclusion 

To summarize, we have proposed that gravitational perturbations about a 
background solution of a unitary and UV-complete theory should be de- 
scribed by an effective theory of a single massless graviton whose linearized 
equations of motion contain at most two derivatives. Given some higher- 
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derivative corrections to Einstein gravity, we have explained how to find this 
theory by imposing a novel set of boundary conditions that ensures uni- 
tarity and removes any unwanted dependence on the cutoff-scale physics. 
In particular, our procedure assigns an unambiguous value to rj/s for any 
higher-derivative model. 

For the purpose of calculating rj/s, there are only two possible candidates 
for a two-derivative description of higher-derivative theories; either Einstein 
or Gauss-Bonnet gravity. We have provided an explicit criterion for choosing 
among the pair that is based on the polarization dependence of their 1PI on- 
shell amplitudes. It has been verified that the resulting value of rj/s agrees 
with our earlier prescription and with previous calculations in the literature. 
We have stressed the issues with directly calculating quantities like r) and s 
in a truly higher-derivative framework. 

To determine the 4- and higher-point functions, in general, further in- 
formation besides polarization dependence versus independence is needed. 
However, if one is satisfied with their leading-order values, then the leading- 
order correction to rj/s is sufficient to determine all higher-point functions to 
leading order, as explained in the last paragraph of Section 3. 

We anticipate that our proposals can be compared with properties of 
strongly coupled fluids such as the quark-gluon plasma or cold atoms. As 
discussed in [29J EO], the LHC should be able to measure energy correla- 
tion functions for a strongly coupled fluid. If this fluid has a gravitational 
dual, these correlators amount to graviton scattering amplitudes. Then the 
polarization sensitivity of the Gauss-Bonnet three-point function should be 
manifest through an angular dependence that is absent in the Einstein case. 
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A The Gauss-Bonnet Lagrangian 

In this appendix, we verify Eqs. ( )25|) and ( )26|) in the main text. 

Let us begin by considering just the four-derivative part of the (general- 
ized) Gauss-Bonnet Lagrangian in Eq. (|23|) . 

pc AGB = pl 2 [n 2 - 4n ab n ab + n abcd n abcd ] , (32) 

where /3 = /3(r) is unspecified and we have dropped e. 

Inserting this Lagrangian into the x, y component of the generic field 
equation (jl]), we have 

V a V b /3 = -g X y/3£ A GB ■ (33) 
This form follows from a special property of any Lagrangian of the 
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Lovelock class (for a formal proof, see Appendix 3 in 

^1 = • < 34 > 

.0 ><-abcd_ 

The next step is to linearize Eq. (1331) and then extract the two- derivative 
terms. The linearization of the right-hand side produces a term proportional 
to h xy , so that this is just a mass term. What remains simplifies greatly in 
the transverse-traceless gauge, as is appropriate for the h xy gravitons once 
the radial gauge has been imposed. The left-most term can be evaluated in 
a tedious but straightforward manner. It can then be expressed as 



pgnh 



xy ; 



where we have defined 

g = 4L 2 



D-2 



K K x x (I) - A)K xy t 



(35) 



(36) 



2D ' '~ 

For sake of completeness, the above expression follows from the kinetic 
part of the linearized field equation or Q ah \3h a b , for which 



( 1 {<c,y} {x,y} \ {x,y} 

c cj^d I c 



be 



(37) 



To obtain Eq. ( 13 6 p . we have set 2Q = Q xy + Q yx and then used symmetries 
of the background to simplify the form of Q. 

Meanwhile, the second term in Eq. (l33|) works out to be 

2eL 2 [g ah Uh xy - V a V b h xy ] V a V b (3 . (38) 

We now want to "reverse engineer" so as to produce the desired result in 
Eq. (|2ip. So let us fix 

G(3Uh xy + 2L 2 [g ab Uh xy - V a V b h xy ] V a V 6 /3 = bU E h xy . (39) 
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To solve the previous equation for the unknown function 0(r) , one can 
impose the hydrodynamic limit, whereby only the radial differentiations of 
h xy turn out to be relevant. This is enough to eliminate the VV/3 term. We 
also take the near-horizon limit, so that □#—>•□. This leads to 

gpnh xy = lnh xy (40) 

or 

(3 = g-'b . (41) 
Finally, inserting the horizon value of —2(D — 1)(D — 4) for Q, we arrive 

at 

g = - 2 (J -ll( C -4) 5 - (42) 

as claimed. 



B Review of the scaling properties of 77 and s 

Here, we review an argument from [28] that is referred to in the main text. 

Let us consider some generic theory of gravity X in a .D- dimensional AdS 
spacetime with a black brane solution. Our interest is in the entropy density 
Sx, which can be obtained by way of Wald's formalism [9, 10J. Then, via the 
analysis of [26J , sx can be expressed in the form of the area law, 

" = ifc • («) 



where a ~ rj^ 2 is the area density of the brane |^| and Gx is related to 



9 We have left a unlabeled, as it is assumed to be fixed for the purpose of comparing 
theories. This is equivalent to comparing theories at fixed or fixed temperature. 
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but generally different than the D-dimensional Newton's constant Ge (the 
subscript E labels quantities in Einstein's theory). 

Now, as shown in [28], there exists a constant conformal transformation 
of the metric (along with an accompanying rescaling of Gx that is necessary 
to preserve the form of the Einstein term in the action) such that 

s ^ 4 • (44) 

with X denoting the transformed theory. That is, the conformal transforma- 
tion is chosen specifically to ensure 

G x -> Gx = G E ■ (45) 

That the entropy density can be freely changed may seem odd, but this is 
because densities such ambiguously defined. Conversely, the entropy 

S = sV± does not change under such a transformation (V± is the transverse 
volume of the brane). Similarly, the ratio of two densities cannot change, 
and so 

2* = ^ . (46) 

The advantage of using the transformed system is that now = se (in 
units of fixed temperature), and so the net correction to the ratio f]xl s x ls 
entirely captured by the viscosity. 

C Review of the Lovelock calculation 



Finally, we recall from [16] why higher-order Lovelock terms cannot con- 
tribute to the ratio r]/s. 
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Let us first consider the entropy density. It is already known that the 
m th -order term C m in the Lovelock Lagrangian (where Cq is a constant, C\ 
is Einstein, £2 is Gauss-Bonnet, etc.), makes a contribution to the entropy 
or its density that, up to constant factors, goes as £ m _i[g||] [Ml [35]. Here, g\\ 
has a special meaning: One should incorporate only curvature components 
of the form TZ XlX ^ iXj (where Xj and Xj are any orthogonal pair of transverse 
brane directions) and permutations thereof. Notably, any such curvature 
component vanishes at the horizon of a black brane; so that the only non- 
vanishing contribution to s is from the m = 1 Einstein term, for which 
An-itsil] = is a constant. 

By direct analogy, the contribution to 77 from the order-m term goes as 
£>m-i\g±] [IS], where g± indicates that one should now only incorporate the 
curvature component lZ rt rt . But, since the background geometry only varies 
with r, one is effectively calculating the m — 1 Lovelock term for a two- 
dimensional theory. As is well known, the p th -order term in the Lovelock 
Lagrangian vanishes identically when p > D/2 . Then, since the dimension 
is effectively D = 2, we have no contribution to rj when m — 1 > 2/2 or 
m > 2 . Meaning that, for the shear viscosity, only the Einstein and Gauss- 
Bonnet terms are contributors. (See, e.g., [36J for an explicit verification.) 
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